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With nonrelativistic QCD factorization, we present the first complete next-to-leading order study 
on the polarization of prompt J/ip hadroproduction by including feeddown from Xc^Pji^ Si) and 
tp{2s) which turn out to be very important parts. By using the color-octet long-distance matrix 
elements obtained from a combined fit of the measurements at the Tevatron and LHC for J/ip, V'(2s) 
and Xc, the prompt J/tp polarization predictions are presented, and the results are in agreement 
with the CDF Run I data (except two points), but in conflict with the CDF run II data, while they 
are close to the ALICE data (inclusive J/'ip). The measurements at the LHC are expected to clarify 
the situation. 

PACS numbers: 12.38.Bx, 13.60.Le, 13.88.-|-e, 14.40.Pq 



Quantum chromodynamics (QCD) is a successful the- 
ory to describe strong interaction, but its fundamen- 
tal ingredients, the quarks and gluons, are not observed 
freely and must hadronize eventually. The fact makes it 
impossible to calculate any processes involving detected 
hadrons in the final or initial states directly. Therefore a 
suitable factorization scheme to divide problems into per- 
turbative calculable and nonperturbative parts is very 
important. In order to remedy the large discrepancy 
between theoretical prediction and experimental data 
on the transverse momentum (pt) distribution of J/ip 
production at the Tevatron, a color-octet (CO) mech- 
anism [l[ was proposed based on nonrelativistic QCD 
(NRQCD) factorization 0, which is proven in next-to- 
next-to-leading order analysis 'i^. It allows consistent 
theoretical prediction to be made and improved pertur- 
batively. However, the leading-order calculation based 
on NRQCD predicts a sizable transverse polarization for 
J ftp hadroproduction at high pt [Jj while the measure- 
ment at the Tevatron Q gives almost an unpolarized re- 
sult. Many efforts have been made to clarify the situa- 
tion thereafter (see Ref. @ and references therein). It is 
expected that the long-standing J/ip polarization puzzle 
could be solved or clarified by the experimental measure- 
ments at the LHC and progress in theoretical calculation. 

In the last five years, there was some very important 
progress made in the next-to-leading order (NLO) QCD 
correction calculation. The NLO corrections to color- 
singlet J7'0 hadroproduction have been investigated in 
Refs. and its pt distribution is found to be enhanced 
by 2 — 3 order of magnitude at the high pt region. And 
it is found that J/ip polarization changes from trans- 
verse into longitudinal at NLO [8]. The results are re- 
produced at leading-order in a new factorization scheme 
for large pt quarkonium production Q . The NLO correc- 
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tions to J ftp production via 5*- wave CO states (S'o^Sp) 
are studied in Ref. [13] and the corrections to pt distri- 
butions of both J/tP yield and polarization are small. 
In Refs. [ll|, NLO corrections for Xc.j hadroproduction 
are studied. The complete NLO calculation for prompt 

3 \S\ 

J I ip hadroproduction (with Pj included) is given by 



two groups [ij, , and their predictions for pt distribu- 
tions agree with the experimental measurements at the 
Tevatron and LHC. Similar progress is also achieved for 
T hadroproduction^^ and J/^p photoproduction Ts"]. 

Recently, the complete NLO calculation for polariza- 
tion of direct J/-0 ( Pj ' included) hadroproduction are 
presented by two groups Jl7j . One of them [iG*! find 
that the CDF data (run II) Q cannot be interpreted by 
using a global fit of J/^ yield data, while the other [13] 
finds that it can be interpreted by using a combined fit of 
J/tp yield and polarization data from the CDF. However, 
there only exist polarization measurements for pro mpt 
(or even inclusive) J/tp production until now [sl. Il8|. It 
is known that among all the feeddown contributions to 
prompt J/^p production from higher charmonium states, 
Xcj contributes more than 20 — 30% of prompt experi- 
mental measurements [H,!!^], and tp{2s) also contributes 
a small fraction, while others are negligible. The feed- 
down contribution is so large that it can drastically 
change the polarization results and must be considered. 
Therefore, to test NRQCD factorization and solve (or 
clarify) the long-standing J/^p polarization puzzle, it is 
a very important step to achieve the polarization predic- 
tions for prompt J/ip hadroproduction. 

In this Letter, we present the calculation of the polar- 
ization for prompt J/ip hadroproduction at NLO QCD. 
In addition to calculating the polarization of J/^ from 

the same intermediate states (^s\^\ ^sf\ ^f^) as 

m previous works we have also calculated 

that from Xcj Cpf\ feeddown, in which the cal- 

culation of the '^j^' part is even more complicated than 

3 fSl 

that of the Pj part of J/tp. In the presented results, the 
polarization for tp(2s),Xcj and prompt J/ip are obtained 
for the first time. 

In NRQCD formalism, the cross section of h hadropro- 
duction is 

(t[pp hx] X! / dxidx2Gl,Gj,a[ij (cc)„a;](0^),(l) 
where p is either a proton or antiproton, the indices i,j 
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run over all the partonic species and n represents the cc 
intermediate states Cs^i \ ^sf\ ^s]^\ ^Pf^) for J/^jj and 

tp{2s), or {Pj , S[ ) for Xcj- The short-distance contri- 
bution a can be perturbatively calculated and the long- 
distance matrix elements (LDMEs) (O^) are governed by 
nonperturbative QCD effects. 

The polarization of J/ ip is described in Ref. [2lj as 



^ _ dgu - do-Qo 
dan + do-Qo ' 



\/2Rcdcr 
ddii + dao 
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where das^st is the spin density matrix of J/ijj hadropro- 
duction. For experimental measurement, a simultaneous 
study of the polarization variables in different reference 
frames is particularly interesting since consistency checks 
on the results can be performed, and it provides a much 
better control of the systematic effects due to detector 
limitations and analysis biases. Therefore, both the he- 
licity and Collins-Soper (CS) frames are employed in our 
polarization calculation. To obtain dcrs^s' i similar treat- 
ment as in Ref. [l5[ is taken for direct J/ip and ip{2s) 
production. The feeddown from 'ip{2s) is obtained from 



da 



J/i' 



S^S'li'i^s) = dagls/B{ij{2s) ~> J/i)). XcO, Xd, 



and 



Xc2 have different masses and decay branching ratios to 
J/'4>, thus are treated differently from J/^{Py). The 
feeddown can be expressed as 



{J/ijiS,hiL,)\xMJz))=ajC: 



(2) 



S L 

Here, Cj j " is the Clebsch-Gordan coefficient, and aj 
is supposed to be independent of Sz, L^, and as an 
approximation. Then the branching ratio of Xcj to Jf-ip 
is approximately expressed as 



E 



(3) 



from which we obtain |ajp. Thus, the spin density ma- 
trix of J /ip from Xcj feeddown is obtained as 



da 



X C 



J, J;: 



J,, J I 



(4) 



The newly upgraded FDC package 122*] is used in our 
calculation, in which the_reduction method for loop inte- 
grals proposed in Ref. [23j is implemented. 

In our numerical calculation. The parton distribution 
function CTEQ6M ^24] and the corresponding two-loop 
QCD coupling constant Us are used. The charm-quark 
mass is chosen as rric = 1.5 GeV and an approximation 
Mh = 2mc is made to fix the masses of quarkonia. The 
color-singlet LDMEs are estimated by using a potential 
model result i25|], which gives \Rj/^{0)\^ = 0.810 GeV^, 

I^V'(2^)(0)P = 0.529 GeV^ and \R'^^iO)\'^ = 0.075 GeV^ 
respectively. Branching ratios are B{J/ip[ilj{2s)] — >■ 
fifi) = 0.0593(0.0077), B[^{2s) J/tp] = 0.595 and 



BiXcj JNA = 0.0116,0.344,0.195 for J = 0,1,2, re- 
spectively [2i|. The factorization, renormalization and 
NRQCD scales are chosen as fir — fJ- f — \/ 4™^ + Vt ^-^id 
/iA — fTT-c, respectively. The center-of-mass energies are 
1.96 and 7 TeV for the Tevatron and LHC, respectively. 
It is well known that the uncertainties for pt distribution 
of charmonium hadroproduction from rric, fJ-A, fJ-r and ^/ 
are large at small pt region. A recent work on relativistic 
corrections to J/ip hadroproduction also shows that 
the correction is negative and large when pt < 10 GeV. 
It is also very clearly shown in Refs. 10, 12] that the ex- 
perimental data in small pt region can not be interpreted 
well. In Ref. [l^, the data in small pt region are also 
included in their fit, but the experimental data in large 
Pt region are sacrificed. Therefore, data in < 7 GeV 
region are excluded in our fit. 

In Ref. [l3| , a combined fit of the prompt J/ip yield and 
polarization measurements of CDF (run II) is performed. 
Without the feeddown contribution, we can repeat their 
results by an exactly same fit. With the feeddown con- 
tribution, we perform a similar fit and find that the pt 
distribution of transverse polarized production rate for 
direct J/ip part becomes negative when pt > 24 GeV, by 
varying the values of fitted CO LDMEs within the un- 
certainty range, the negative behavior of direct part can 
be delayed as pt increases, but cannot be avoided, which 
means that no physical solution in the LDME parameter 
space can be achieved. It clearly shows that the feed- 
down contribution is so important that the conclusions 
for the polarization arc completely different for the cases 
with or without the feeddown part. Therefore, we choose 
to fit the CDF and LHCb experimental data for the yield 
only and give predictions for polarization in this work. In 

addition, p'^^^ ~ p^{Mj/^/Mh) is used as an approxima- 
tion to count for the kinematics effect in the feeddown of 
h = iIj{2s),XcJ- 

For XcJ production, the CO LDME is obtained 
(O^-^fsf^)) = (0.221 ± 0.012) X 10-2 GeV^ with 
2.57/10 by fitting the CDF and LHCb 
data [ig|,l20|,[2^. 

The CO LDMEs of tp{2s) are obtained from a com- 
bined fit of the CDF and LHCb experimental data [131 



as 

xVd.oi 



((o(M'')),(o(M")>, 



{o(pr)>s 



O = (-0.012±0.869,0.34±0.12,0.42±0.24) 

with x^/d.o.f. = 17.8/17. The large uncertainty arises 
from approximately linear dependence of three short- 
distance parts. To clarify the situation, a rotation matrix 
V as discussed in Ref. [12] is used to make fitting variables 
independent by introducing A = OV, and A is obtained 
with only independent error for each A^ in the fit. Then 
differential cross section da is obtained with 

da = Y^ 0,da, = ^ OVV-^da = ^ AV-^da. 

A = (0.16±0.91,0.47±0.05,-0.21±0.01), (6) 
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FIG. 1: pt distribution of prompt J/^p and "0(28) hadroproduction. CDF and LHCb data are taken from Refs. [2^ - l30l ]. 
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FIG. 2: Polarization parameter A of ■(/'(2s) in helicity(up) and CS(down) frames. The CDF data are taken from Ref. [sl. IT^ 



where A is obtained with same x^/d.o.f. but much 
smaller uncertainty. With the values of A and relation in 
Eq. ([S]), theoretical predictions with proper uncertainty 
can be made. 

After the treatment of feeddown, the CO LDMEs for 
J/'0 are obtained from a combined fit of the CDF and 
LHCb experimental data [1^ [2^. By using the same 
definition in Eq. ([5]) for J/ip, the fit gives 

O = (9. 7 ±0.9, -0.46 ±0.13, -0.95 ±0.25), 

A = (-9.6± 1.0,1.7±0.1,-0.37±0.01), (7) 

with x^/d.o.f. = 5.32/10. Thereafter, we use the error 
of each independent variable in A to generate all the un- 
certainty bands in theoretical predictions. It is clearly 
shown in the following figures that the uncertainty band 
is not too wide even with large uncertainty in A in Eq. ([5]) . 

In Fig. [1] we find that the feeddown part contributes 
almost the same as the direct part in prompt J/^p yield 
when Pt > 25 GeV at the Tevatron, so as for \y\ < 2.4 
when Pt > 30 GeV at the LHC, but is less important in 
the forward range 4.5 > y > 2. The V'(2s) polarizations, 
shown in Fig. [21 go from longitudinal to transverse as 
Pt increases in the helicity frame, which has totally op- 
posite trend with current CDF measurement 5], and go 



from transverse to slight longitudinal in CS frame. As 
is shown in Fig. |31 the polarizations of J/ip from Xcj 
feeddown result in small transverse polarization 0.2) 
in the helicity frame, and slight longitudinal polarization 
—0.1) in the CS frame at large pt- Finally our the- 
oretical predictions for prompt J/ip polarization in both 
helicity and CS frames are shown in Fig. 31 in comparison 
with current existing measurements from the CDF and 
ALICE Collaborations. In the forward rapidity region, 
our predictions are close to the ALICE measurement for 
inclusive J/ip production in both helicity and CS frames. 
In the central rapidity region, our results are in agree- 
ment with the CDF run I data (except two points), but 
in conflict with the CDF run II data. However, we still 
cannot draw a definite conclusion since there is no way to 
judge these two measurements. From this point of view, 
this is another reason for us to exclude CDF data on J/ i/j 
polarization in our fitting. With the optimized analysis 
method to measure more J/^ polarization information 
with two different frames (helicity and CS) as used in 
the ALICE measurement, it is expected that J/ip and 
tpi'is) polarization measurements at the LHC would help 
to solve the polarization puzzle or clarify the situation. 

The polarization predictions in Ref. [l^ and [l3| are 
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FIG. 3: Polarization parameter A of J/i/) from Xcj feeddown in helicity(up) and CS(down) frames. 



for direct J/?/; production. Comparing them with the 
CDF (ALICE) measurements for prompt (inclusive) J 
production and making a definite conclusion can be con- 
sidered as reasonable when and only when the feeddown 
contribution is negligible. Our calculation shows that the 
feeddown contribution is very important and cannot be 
neglected. Our predictions are in better agreement with 
the ALICE measurement than theirs in Ref [l^. By ex- 
cluding the feeddown contribution, our polarization pre- 
dictions for direct J/-)/) production are consistent with the 
results in Ref. |16] and [17] by using their fitted LDMEs. 

In summary, we presented the first complete NLO 
study on polarization for prompt Jj^l) and ipi'^s) 
hadroproduction in both the helicity and CS frames by 
including feeddown contributions from higher charmonia 
Xc and ?A(2s) based on the NRQCD factorization scheme. 
It is found that there is no physical solution for a com- 
bined fit of both J ftp yield and polarization measure- 
ments by CDF at the Tevatron (run II) for prompt J/V' 
hadroproduction, while the previous fit without con- 
sidering the polarization from feeddown of Xc and 'ip{2s) 



can be reproduced. It means that the polarization con- 
tributed from Xc and ip{2s) feeddown is a very impor- 
tant part and must be included to solve or clarify the 
J polarization puzzle. Therefore, we choose to use the 
new CO LDMEs from a combined fit of J/ip yield mea- 
surements at the Tevatron and LHC with pt > 7 GeV. 
And the NLO theoretical predictions on polarization for 
prompt J/ip are presented at the Tevatron and LHC. The 
results are in agreement with the CDF run I data (except 
two points), but in confiict with the CDF run II data. It 
is close to the ALICE measurement at the LHC although 
the measurement is for inclusive J/ip. Finally, it is clear 
that the polarization measurements at the LHC are very 
important to clarify the long-standing J/^p polarization 
puzzle. 
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